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We consider a Boussinesq system describing one-dimensional internal waves 
which develop at the boundary between two immiscible fluids, and we restrict 
to its traveling waves. The method which yields explicitly all the elliptic or 
degenerate elliptic solutions of a given nonlinear, any order algebraic ordinary 
differential equation is briefly recalled. We then apply it to the fluid system 
and, restricting in this preliminary report to the generic situation, we obtain 
all the solutions in that class, including several new solutions. 
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1. Introduction 

At the boundary between two immiscible fluids, one observes the formation 
of waves, called internal waves. These are typically described by Boussinesq 



* Waves and stability in continuous media, eds. A. Greco, S. Rionero and T. Ruggeri 
(World scientific, Singapore, 2010). WASCOM 15, Mondello (Pa), 28 June-1 July 2009. 
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systems such as 2 

{(1 — [ibA)d t ( + c v div v + cnl div(£v) + jj,a div Av = 0, 
(1 — ndA)d t v + c z grad£ + -^p grad |v| 2 + itcAgrad( = 



(1) 



in which /i, a, &, c, d, cnl, c v , c z are constant, /icnl ^ 0. 

We restrict here to one-dimensional situations, relevant for instance 
when the fluids are inside a channel, and our purpose is to obtain traveling 
waves C = — c v/cnl, v = v(£)+co/cnl,£ = x — cot (the translation on u 
suppresses any dependence on c v and the one on v shortens the expressions 
below) in closed form by a nonperturbative method. The conservative form 
of the equations (1) allows each equation to be integrated once, and the 
considered system will be 



The above system is essentially the same as those considered by Chen 4 
and Nguyen and Dias 10 . 

Chen 4 already found all the traveling waves in which u and v are 
polynomials (of degree 1, 2 or 4) in tanh fc£ and sech k£. In the present work, 
we obtain the closed form expressions of all & those solutions of (2) which 
arc cither elliptic (doubly periodic in the £ complex plane) or degenerate 
elliptic, i.e. rational in one exponential e k ^ (simply periodic in the £ complex 
plane, which includes the above mentioned solutions) or rational in £. 

The method, based on those complex singularities of (u, v) which depend 
on the initial conditions ("movable" singularities 5 ), implements classical 
results by Briot, Bouquet 3 and Poincare 7 . First presented in 9 , it was later 
turned into an algorithm 6 . 

2. Singularity analysis 

In order to know whether closed-form solutions to (2) might exist, a pre- 
requisite 5 is to investigate the singularities of (u, v) in the complex plane 
of £. One must distinguish whether the total differential order of system (2) 
is four or two, depending on the value of the determinant 



We leave the nongeneric case S — to a forthcoming detailed study. 





S = bdc\ — ac. 



(3) 



a In this short report, only the generic case is presented. 
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In the generic case 5 ^ 0, the system is equivalent to 

^Su" + c NL (dc uv — t^ 2 ) — ac zU + dc a Ki — aK 2 = 0, 



Ll5v" + cnl ( ^% 2 — cuv ) + bc c z ii + bc$K 2 — cK\ = 0. 



(4) 



For an easier computation of u knowing v, it is convenient to introduce the 
shift 

bc ° 

v s = v c z . (5) 

cnlc 

In order to find all the elliptic and degenerate elliptic solutions, one must 
first determine all the families of movable poles and movable zeros (i.e. mov- 
able poles of 1/u or l/v s ). 

Let us first determine the poles. Assume that, near a movable singularity 
£o, the variables (u, v) behave algebraically 

u ~ u x Pl , v ~ v X P2 , X = £ - Co, uovo ^ 0, (6) 

with pi,^2 not both positive integers. Balancing the highest derivatives and 
the nonlinear terms, one generically (nongeneric cases will be dealt with in 
a forthcoming paper) obtains double pole behaviours 

f Pl = -2, p 2 = -2, 

| c NL u v + 6fi(bc u + av ) = 0, -4p«o + 6v-(dc v + cu ) = 0. 
Whenever acD ^ 0, with 

D 2 = {b- 2d) 2 c 2 + Sac, (8) 
this system (7) admits two solutions (u ,v ), 



(7) 



SacD ^ : < 



UQ.e = 7T^— ( b ( 2d - b ) C l - 4flC + £bc D ) > 

2cc NL 

vo.e = — (-(& + 2d)c Q + eD) , e 2 = 1. 
cnl 



(9) 



Let us next determine the movable zeros. By elimination between (4), it is 
easy to establish the fourth order ODE for l/v s . Its movable poles are one 
simple pole of arbitrary residue plus, when K\\ ^ 0, one double pole, with 

Ku =Kl - XD - c 2dc ° [k 2 + ^(XD- 2dc ) 2 ) . (10) 

The movable zeros of u are less easy to establish, but it is sufficient for 
our purpose, as explained below, to know that u has always at least two 
movable zeros. 
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One must then compute 5 the Fuchs indices i of the linearized system of 
(4) near all the movable singularities. Near the movable double poles (6), 
the resulting indicial equation (we skip the details) only depends on one 
adimensional parameter A, 

(* + 1)(» - 6) {f-^ + YT^xj =0 ' 

A= (6+ ^ )C0 . (12) 

Since some Fuchs indices i are generically noninteger, the general solution 
of the system (4) is multivalued. Nongenerically, for the general solution to 
be singlevalued, it is necessary that, for both signs e, all roots i of (11) be 
integer. Denoting these roots as 

1 . 5 + Q 5-Q 5 + i? 5-R 

- 1, 6, — ^— , — — for e = +1, -1, 6, — — , for e = -1, (13) 

the elimination of A between the products of the roots 
25 - 2 12 25 - R 2 12 



4 1 + A' 4 1-A' 

yields the diophantine equation 



(14) 



24 24 
25~ 1 + 25"^ = lj (15) 

which admits no solution for odd positive integers (0, R)- 

Despite its generically multivalued general solution, the system (4) may 
still admit singlevalued particular solutions. For this it is necessary that the 
Laurent series whose first term is (6), 

u = ^2 u jX j ~ 2 , v = ^v j x 3 ~ 2 , (16) 

3=0 3=0 

exists, i.e. that no impossibility occurs when computing the coefficients 
Uj , vj . The invariance of the system (2) under £ — > — ^ forbids the occurrence 
of odd powers of \ m the Laurent series of u and v. 

For convenience, the three defining equations for S, D, A can be solved 
for a, bca, S in terms of dco, A, D, c, yielding 

D 2 - (Adc - \D) 2 



8c 



bcO = -2dc + AD, 5= (A - l)D 2 , (17) 



u , e = P^-(X - e) [Mc - D(X - e)] , v , e = ^-(e- A). (18) 

4CCNL CNL 
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For a generic A, the values i = 2 and i = 4 arc not roots of (11), 

( b + 2d ^o 4 {1 4} (19) 
(b - 2dfc 2 + Sac * 1 ' j ' U j 

so no impossibility can occur when computing the next coefficients 
U2,v 2 ,U4,V4. The only obstruction arises from the Fuchs index i = 6, which 
generates two necessary conditions (one for each sign e) for the absence of 
movable logarithms (again we skip the details of this classical computation), 

Q 6 = c z {XD - 3dc ) (kx-^^-k^] -0, (20) 
V cv J 

2 2 

K x = cKx - a(b - 2d) C (A, K 2 = K 2 + ({b - 2d) 2 c 2 + 2ac) (21) 
This defines three subcases, 

c z = 0, (22) 

(b - d)c Q = 0, (23) 

Kx = 2 -^K 2 . (24) 

The first two subcases are independent of the sign e, and the third 
subcase can be enforced either for one sign (condition (24)) or for both 
signs, leading to the stronger condition 

Kx - 0, K 2 = 0. (25) 
For the first two cases, a first integral exists, 
(b — d)c c z = : K 6 = yUC NL (c 2 u' 2 + 2cdc u'v' + (ac — (b — d)dc^ ) )v'' 2 ^j 
+ c nl (cuv 2 + fy c ° v 3^j _j_ CNLCCzU 2 

+2cK 2 u + (2cKx - 2K 2 (b - d)c ) v, (26) 

and the first integral for the third subcase kx = K 2 = 0, yet to be found, 
is not quartic in (u' , v'). 



3. Method to find all the elliptic solutions 

For full details on the method, we refer to 6 ' 5 . 

The input is an iV-th order (N > 2) any degree autonomous algebraic 
ordinary differential equation (ODE) admitting a Laurent series. 

The output is made of all its elliptic or degenerate elliptic solutions in 
closed form. 
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Let us first recall a classical definition. The elliptic order of a nondc- 
generate elliptic (genus one) function 1 is the number of poles, counting 
multiplicity of course, inside a period parallelogram. It is equal to the num- 
ber of zeros. This equality breaks down under degeneracy to genus zero, 
e.g. for the rational function u = (£ — a)(£ — — c). 

The successive steps of the algorithm are 6 : 

(1) Find the analytic structure of singularities (in our case two families of 
movable double poles for both u and v, see (7), one movable simple zero 
and, if K\\ ^ 0, one movable double zero for v, at least two movable 
simple zeros for u). Deduce the total number of poles (or, if greater, of 
zeros) of the unknown function and its derivative, here m = 4, n = 6 
for (v,v'), more for (u,u'). 

(2) Compute slightly more than (m + l) 2 terms in each Laurent series. 

(3) Choose one of the dependent variables (u, v) (call it U) and define the 
first order m-th degree subequation F(U, U') = (it contains at most 
(m+ l) 2 coefficients a,j.k), 

m 2m — 2k 

F(U, U') = J2 E a J,kU j U' k - 0, ao.r n ? 0. (27) 

fc=0 j=0 

(4) Require at least one Laurent series of U to obey F(U, U') = 0, 

F ee F oX 3 + 0(x J+1 ) j , Vj : Fj = 0, (28) 

and solve this linear overdetermined system for a^fc. 

(5) Integrate each resulting first order ODE F(U, U') = 0. 

The key advantage of this method is that the system of equations Fj = 
for the unknown coefficients cij.k is linear and infinitely overdetermined, 
therefore quite easy to solve. 

4. Elliptic and degenerate elliptic solutions, generic case 

By generic, we mean that the fixed constants a, c, bco, dco, c z , K\, K<i obey 
the nonvanishing conditions acSD ^ 0, (19) and only one of the three 
vanishing conditions (22), (23), (24). 

When the algorithm of section 3 is applied to a system of ODEs such 
as (2), a key practical ingredient is to select a "good" variable U, i.e. one 
whose total number of poles (or, if greater, of zeros) of U and U' is the 
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smallest possible. Since this number is always smaller for (v,v') ((4,6)) 
than for (u, u'), the natural choice is U — v. 

Moreover, the already mentioned invariance under £ — >• — £ forbids the 
occurrence of odd powers of V = v' in the subequation (27). 

Since U = v admits two distinct Laurent series, the search for elliptic or 
degenerate elliptic solutions splits into (step 4): either require one Laurent 
series to obey (the odd-parity terms have been removed) 

F = U' 2 + a 3 . U 3 + a 2fl U 2 + a lfi U + a ,o - 0, (29) 

or require both Laurent series to obey 

F = U' 4 + a 3 , 2 U 3 XJ' 2 + a 6fi U 6 + a 2 , 2 U 2 U' 2 + a bfi U 5 + a h2 UU' 2 + a Afi U A 
+a oa U' 2 + a 3fi U 3 + a 2fl U 2 + a lfl U + a , = 0. (30) 



4.1. Elliptic and degenerate elliptic solutions, one series 

One of the three necessary conditions (22), (23), (24) must hold true. In 
this section we simply denote (u ,vo) anyone of the two values (uo i£ ,Vo, e ) 
(Eqs. (9) or (18)). 

In step (2) of section 3, it is quicker to compute simultaneously both 
Laurent coefficients (uj , Vj ) from system (4) . Going to j = 8 is sufficient 
to obtain all the coefficients in (29) and ensure that system (4) is indeed a 
differential consequence of (29). 

In step (4), with the definition (5), the resulting subequation (29) is 

«f - — v 3 + b 2 v 2 s + hv s + b = 0, (31) 

cJb - d)c = : &i = 0, b = 0, K x = 2— i^, (33) 

cv 

c z (b - d)c ^ : &i = — !— ( n , 8 .^ 2 , - 4(3dc - \D)% ) , b = arb, 
c N LM \D{\ - e) c 2 J 

i.e. one additional constraint is found in the case c z (b — d)co = 0. 

Step (5) is immediate. Indeed, subequation (31) is nothing else, up to 
an affine transformation, than the canonical equation of Weierstrass 1 , 

v s = v p(£ - £0,32,33) +v 2 - -^-c z , (34) 

cnlc 

p' 2 = 4p 3 - g 2 p - 33 = 4(p - ei)(p - e 2 )(p - e 3 ), p" = 6p 2 - y . (35) 
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The result is 

Du c z v 6 
v = v p- — , 53 = 28—, 

u = ^ f Vs - ($dc - D(3X - e))- 



v \ 2cc NL 



3 (36) 

Cz (6 - rf)c = : K X = ((3e - A)D - 4dc )^, 52 = ^f^, 53 = 
c z (b - d)c a ^ : tYi = 2DU ° K 2 , 53 = arbitrary, 

52 = (dw (-^r + (12dco - (3A + • 

The effective expression of p depends on the values of (32,53) in (35), 
according to the identities 1 , 

( doubly periodic ( "cnoidal wave" ) , g 2 — 
p(z 7 g2,g 3 ) = I 3gcoth 2 (y3^z) - 4q, 52 = 12g 2 ,s 3 = Sq 3 , (37) 
I ~ ' 92 = 93 = 0. 

For each of the three subcases (22), (23), (24), one thus obtains two 
closed form solutions (one for each sign e), namely: 
two rational solutions for the condition (22), 

(22) : v s = v = v (H £o)- 2 , u = - £ ) 2 , 

K x = ((3e - X)D - 4dc )^, v 6 = 0, K 6 = 0, (38) 
two solutions rational in one exponential for the condition (23), 

(23 : v = v [t z - — - -5 , r = - tanh , 



«o V 6cc N l/ 



Dc z \ K 2 v 



M CC NL Us 



and two doubly periodic solutions for the condition (24), 
(24) : v = v p + -^(12dc - D(5X - e)), 

4CCNL 
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^ (v s ~(8dc a -D(3\-e))-^ 

V \ ^CCNL 

-32c 2 K 2 + cl(l2dc - £>(3A + g)) 
= \2(\-e) 2 \i 2 D 2 c 2 



28« 6 - 2Du ~ 
03 = , ili = K.2, vq — arbitrary. (40) 

V CV Q 

When only one Laurent series for v is enforced, the method therefore 
yields two solutions of each possible kind (elliptic, rational in one exponen- 
tial, rational). 

Remarks. 

(1) If one is interested in finding only the nondegenerate elliptic (genus one) 
solutions, such as (40), a quicker method to find them all is to combine 
the present method with the conditions on the residues as explained in 

8 

(2) For the rational and trigonometric solutions, the ODE obeyed by u has 
the type (30) (degree four), e.g. for the rational solution u = a(£ — 

£o)- 2 + /3(£-£o) 2 , 

(au 2 - 2u 3 + 8a/3u) ^ - 4 (u 2 - Aa0) 3 = 0. (41) 
This is why u should not be chosen to apply the algorithm. 

4.2. Elliptic and degenerate elliptic solutions, two series 

One of the three necessary conditions (22), (23), (25) must hold true. 

The subequation (30) is assumed to be nonfactorizable (nonzero value 
for the discriminant in U' ). In order to determine all the aj t k in subequation 
(30) and ensure that system (4) is indeed a differential consequence of (30), 
it is necessary and sufficient to compute the series up to j = 14 included, 
i.e. 8 terms in the series for u and 8 terms in the series for v if one uses the 
system (4) to perform the computation. 

For the condition (25), the found subequation (30) is the product of two 
factors like (29), the solution would only represent (40), so we discard it. 

For either condition (22), (23), one finds the unique subequation 

c z (d - b)c = : F = [v' s 2 + T(v s - A)(v 2 s - B)] ' - A(v 2 s - Bf = 0, (42) 

in which the constants (A, B, T, A) take the values 

3(\ 2 -l)D \ 2 -l - 4c NL A 

4cc NL A c^ h (\ 2 - 2) 3nD(\ A - 1) 
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A = ( ^ V, k x = (-2dc + ^f-f) k 2 . (44) 
\3fiD(X 2 - 1)7 V 2(A 2 -2) y 1 ; 

Because of the parity invariance, the general solution of (42) can- 
not involve p' and is the quotient of two second degree polynomi- 
als of p. It can be obtained by brute force with the Maple command 
with(algcurves) ;Weierstrassf orm( . . .) 7 , but the structure of singu- 
larities allows a straightforward integration. Indeed, the a priori solution 

a 2 p 2 + a 1 p + a 

Vs = 7 \7 T> ( 45 

(p-ei)(p-e 2 ) 

must have two double poles with principal parts v . e (£, — £ £ )~ 2 , therefore 
the poles e\, e 2 must be 1 two of the three zeros ej of p', Eq. (35), 

Pfe) = ej, j = 1, 2; ^ & : p(£) - ej ~ ^ (6e 2 - f ) (£ - (46) 
The result, 

^0,1 f c „ 2 52 N \ 1 , «0,2 ^ 2 52 ^ 1 



vs = a 2 + ^ 6ef - ^ — — + ^ 6^ - ^ — (47) 
2 \ 2/p — ei 2 \ 2 / p — e 2 

in which vo,i, «o, 2 are the two values of «o, e , Eq. (9), can be written as 

3/"£> „. , (2p + e 3 + A(ei -e 2 ))(p-e 3 ) 

w s = 3Ae 3 + ei - e 2 - (ei - e 2 )- 



cnl " " ' " (p-ei)(p-e 2 

ei + e 2 = -e 3 



12 M c' v 1 " 3AV(A 2 -2) z V 4 /W 



(48) 



2 



_ c z if 2 / c. V 2K 2 f c z \ 

ff2 6/xc3AV(A 2 -2) + Umc/ ' 33 3A 2 M 2 (A 2 -2) UaW ' 
This represents two solutions since ei — e 2 can take two opposite values. 

This solution is trucly elliptic, except when two of the three roots 
ei, e 2 , e 3 are equal. Finally, the value of u is easily obtained from the second 
equation of system (4). 



5. Conclusion and perspectives 

When the coefficients of the Boussinesq system (1) take generic values as 
defined at the beginning of section 4, we have obtained in closed form all 
the singlevalued solutions which are either elliptic ("cnoidal") or rational 
in one exponential (this includes polynomials of tanh and sech) or rational. 
Naturally, one then must select those solutions which make the internal 
waves real and bounded. 
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When the quadratic nonlincarities in the Boussinesq system (1) are in- 
sufficient and cubic nonlinearities must be incorporated, the resulting sys- 
tem 10 can be handled similarly. The various Boussinesq systems displayed 
in 2 could similarly be processed. 
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